function dydt=o_PPN(t,y,G,P,k)
% y[1]=pos_x;y[2]=pos_y;y[3]=

